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Abstract 

The conventional definition of the running coupling a^g(fi) in quantum chro- 
modynamics is based on a solution to the renormalization group equations 
which treats quarks as either completely massless at a renormalization scale 
fji above their thresholds or infinitely massive at a scale below them. The 
coupling is thus nonanalytic at these thresholds. In this paper we present an 
analytic extension of a^g(fj,) which incorporates the finite-mass quark thresh- 
old effects into the running of the coupling. This is achieved by using a 
commensurate scale relation to connect a^g(fx) to the physical ay scheme 
at specific scales, thus naturally including finite quark masses. The analytic- 
extension inherits the exact analyticity of the ay scheme and matches the 
conventional MS scheme far above and below mass thresholds. Furthermore 
just as in ay scheme, there is no renormalization scale ambiguity, since the 
position of the physical mass thresholds is unambiguous. 
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I. MOTIVATION 



The running coupling in quantum chromodynamics (QCD) in the modified minimal 
subtraction (MS) scheme |l] and other dimensional regularization schemes is traditionally 
constructed by solving the renormalization group equations using perturbative approximants 
to the f3 function which change discontinuously at the quark mass thresholds f|-|J. This is 
equivalent to using effective Lagrangians with a fixed number of massless fermions in each 
energy range between the quark mass thresholds. Thus in the MS scheme, the /3(/z) function 
depends on the number of "massless" quarks N F (n) which is taken as a step function of the 
renormalization scale /i. Matching conditions at threshold require the equivalence of one 
effective theory with n massless flavors to another effective theory with one massive and 
(n — 1) massless quarks. It should be noted that this does not prevent one from including 
quark masses in the MS scheme. However, the quark masses do not enter into the (3 function 
since the running of the coupling is mass independent. 

The one-loop matching conditions [|]-[7| in the MS scheme require the coupling to be 
continuous if the matching is done at the quark masses, although the derivative is discontin- 



uous. In two-loop matching |p|-|T0[| the coupling itself becomes discontinuous if the matching 
is done at the quark masses, but it can be rendered continuous by modifying the MS scheme 
H]. Recently the three-loop matching conditions have been computed [110], which together 



with the four-loop /9-function |TT|, give the possibility to evolve the MS coupling to four loops 



with massless quarks. This gives a reduced dependence on the matching scale, as shown in 



l2]| , but possibly a nonphysical threshold dependence. However, in such a treatment the 
derivatives of the coupling remain discontinuous. The inevitable result of the matching in a 
dimensional regularization scheme is that the running of the MS coupling in the renormal- 
ization scale is nonanalytic - nondifferentiable or even discontinuous - as the quark mass 
thresholds are crossed. Thus there is an intrinsic difficulty in expressing physical, smooth 
observables as an expansion in the MS coupling. It is clearly necessary to restore the finite 
quark mass effects in their entirety in order to restore analyticity. 

Aesthetically, it is unnatural to characterize physical theories in terms of an artificially- 
constructed renormalization scheme such as MS; it is more physical to use an effective charge 
as determined from experiment to define the fundamental coupling fl3f . For example, in 
analogy to quantum electrodynamics, one could choose to define the QCD coupling as the 
coefficient oty(Q) in the static limit of the scattering potential between two heavy quark- 
antiquark test charges: 

V(Q 2 ) = -4nC F ^l (1) 

at the momentum transfer q 2 = t = —Q 2 , where C F = (N^ — l)/(2Nc) = 4/3 is the 
Casimir operator for the fundamental representation in SU(Nc), (with Nq = 3 for QCD). 
Such an effective charge automatically incorporates the quark mass threshold effects in the 
running, and thus it has an analytic (3 function. The ay scheme is particularly well-suited 
to summing the effects of gluon exchange at low momentum transfer, such as in evaluating 



the final-state interaction corrections to heavy quark production [14], or in evaluating the 



hard-scattering matrix elements underlying exclusive processes [|l5j. A physical effective 
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charge has the additional advantage that the Appelquist-Carazzone decoupling theorem [16 
is automatically incorporated. 

In this paper we shall construct an analytic extension of the «ms scheme, which we call 
5ms, by connecting the coupling directly to the analytic and physically-defined ay scheme. 
The necessity for an analytic coupling has been emphasized by Shirkov and his collabora- 
tors [[L7|]. Our definition allows one to use a scheme based on dimensional regularization, 
but which also, in a simple way, treats mass effects properly between the mass thresholds. 
Thus, instead of having the number of effective flavors (Np) change discontinuously at (or 
nearby) the quark threshold, we obtain an analytic Np(fi) which is a continuous function 
of the renormalization scale /i and the quark masses m^. Thus the analytically-extended 
scheme inherits the mass dependence of the physical scheme. In addition, the renormaliza- 
tion scale \i that appears in the analytically-extended scheme 5m§- is directly related to the 
momentum transfer appearing in the ay scheme and thus has a definite and simple physical 
interpret at ionF|. 

The essential advantage of the modified scheme 5ms is that it provides an analytic inter- 
polation of conventional dimensional regularization expressions by utilizing the mass depen- 
dence of the physical ay scheme. In effect, quark thresholds are treated analytically to all 
orders in m 2 /Q 2 ; i.e., the evolution of our analytically extended coupling in the intermediate 
regions reflects the actual mass dependence of a physical effective charge and the analytic 
properties of particle production in a physical process. Just as in Abelian QED, the mass 
dependence of the effective potential and the analytically-extended scheme 5m§- reflects the 
analyticity of the physical thresholds for particle production in the crossed channel. Fur- 
thermore, the definiteness of the dependence in the quark masses automatically constrains 
the renormalization scale. Alternatively, one could connect «ms to another physical charge 
such as an defined from e + e~ annihilation. 

Our approach should be compared with the standard treatment of quark mass threshold 
effects in the MS scheme. For fixed order in a s the corrections due to finite quark mass 
threshold effects which we are considering in this paper have been calculated for the hadronic 



width of the Z-boson and the r lepton semihadronic decay rate [|18H20|J9|1 . The calculations 
have been made both exactly to order a 2 ^ and as expansions in terms of m 2 /Q 2 and Q 2 /m 2 
for light and heavy quarks respectively. Note that in principle the determination of the finite 
mass threshold effects for physical observables in dimensional regularization schemes would 
require a complete all-orders analysis of the higher-twist mass corrections to the effective 
Lagrangian of the theory. 

There are a number of other reasons to construct an analytic extension of the «ms scheme: 

• The comparison of the values of the coupling a s as determined from different experi- 
ments and at different momentum scales is an essential test of QCD (for a recent review 



of existing measurements see |2lfl). One source of error is neglect of quark masses in 
the determination of a s and in the subsequent running of the coupling from the scale 
where it has been determined to the conventional reference scale, the Z-boson mass. 



1 A somewhat similar approach has been tried in [17], but using the unphysical MOM renormal- 
ization scheme to implement the mass thresholds. 
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Lattice calculations for the J/^ and T spectra now provide the most precise deter- 
mination of a s at low momentum scales [22-25]. It is important to know how finite 



quark mass effects enter into the running of this value of a s to lower and higher energy 
scales with as small an error as possible. 

Finite mass threshold effects in supersymmetric grand unified theories are important 
when analyzing the running and unification of couplings over very large ranges. It has 
been discussed, for example, in refs. [p6,|27|. However, the scale used in the running and 



for the threshold effects has not been related to the physical scale which is naturally 
obtained in our approach. 

It is natural to unify theories by matching physical couplings and masses at the unifi- 
cation scale. This can be accomplished in the ay scheme or equivalently 5ms. 



II. DETAILS OF ay 



In the case of the Abelian theory, the coupling ay derived from the heavy lepton potential 
is equivalent to using the effective charge defined from the running of the photon propagator. 
In the non- Abelian theory the gluon propagator is not gauge invariant; one thus has to turn 
to a physical gauge invariant observable such as the heavy quark potential. 

The effective charge ay(Q), defined as in Eq. (jl|) can be calculated as a perturbation 
expansion in a-M§, 



a v (Q) = «Ms(/i) + v lt i 




+ v 2 ,i 




(2) 



The first two nontrivial terms in the perturbative series have been computed in the MS 
scheme [28-33]. A comprehensive analysis of ay to order a 3 has recently been given by M. 
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where (3$ and /3\ are the first two universal coefficients in the /3-function. It is also known 
that the next coefficient in the expansion is non-analytic in a = 0, since it contains a ln(a) 
term |[30|| . This non-analyticity does not originate from fermionic corrections to the heavy 
quark potential. This can be seen by adopting a physical gauge. In such gauges we would 
reproduce the same analyticity structure for fermionic corrections in QCD as we find in 
QED. In QED, however, we have no problem with analyticity at any order in perturbation 
theory. Note that the total derivative of ay with respect to the renormalization scale /x of 
the «ms scheme is zero, since ay is a physical observable. 
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The scale (Q) dependence of the effective charge defines the equivalent of the Gell-Mann 
Low ^-function for the effective charge in QED Q. In the case of ay, Q is the momentum 
transfer in the heavy quark potential, and the -^-function is given by, 



da v (Q) 
dlnQ 
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7T 
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(4) 



The first two terms in this series coincide with the universal and well-known first two terms 
in the Callan-Symanzik /^-function [35|, i.e. ip^ = (3q and ip^ 1 ' = Pi, but the higher order 
terms, ipy etc., depends on the observable under studyf]. 

For completeness we give the coefficients in the QCD ip function with the normalization 
used above, 
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The results given above in the scheme have been obtained using massless QCD. The 
effects of non-vanishing quark masses can be taken into account by using a (^-dependent 
N F which will be derived from the one-loop massive vacuum polarization function in the 
next section. 



III. CALCULATION OF THE RUNNING COUPLING TO ONE-LOOP ORDER 

Our approach in this paper is as follows: the ay{Q) scheme automatically includes the 
effects of finite quark masses in the same manner that lepton masses appear in Abelian 
QED. We can then relate the MS scheme to the ay scheme through a commensurate scale 



relation p6f, which is effectively a scale transformation between the two schemes. The 



2 Some authors denote the coefficient ifiy by Pzy instead. The convention used here is to empha- 
size the difference between the dependence on the physical scale Q and the unphysical renormal- 
ization scale /i. 
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analytic dependence of ay is then transferred to the analytically-extended 5ms scheme. The 
usual massless expressions are recovered far above or far below any individual quark mass 
threshold. 



A. Calculation of the Mass Dependence for the Running Coupling 

The coupling ay(Q), which is derived from heavy quark scattering, is closely related to 
the renormalization of the gluon propagator. In physical gauges with Z\ = Z2 the coupling 
renormalization is due purely to self-energy insertions in the propagator.^ 

For the purposes of this paper it will be sufficient to restrict our analysis to one-loop 
ordeif], i.e. 

The physical running coupling in the ay scheme, normalized at an arbitrary momentum 
transfer scale Qo, may be represented as 

a v (Q) = . (5) 

l-U(Q,Q ,a v (Q )) 

The vacuum polarization function IT may be computed from the perturbative expansion of 
the renormalized propagator between heavy quarks. The coupling is then 

av(Q) = a v (Qo) [1 + n + n 2 + n 3 + . . .] , (6) 

where we have used the shorthand II = H(Q, Qo, ay{Qo)) for the renormalized sum of all 
one-particle irreducible 1PI diagrams for the gluon self-energy. Since the coupling has the 
value ay{Qo) = a at the physical renormalization point Q = Qo, the self-energy obeys the 
boundary condition Tl(Qo, Qo, Q^o) = 0. 

We begin by considering the integral representation of the quark part of the one-loop 
gluon vacuum polarization diagram (see Fig. p]): 



3 Strictly speaking, this is only true up to one-loop in QCD and two-loops in QED. At higher 
orders new types of diagrams appear in the potential which cannot be described as simple self- 
energy insertions in the propagator. In QCD such a diagram is the so called "H-graph" and 
in QED the light-by-light scattering diagram has the same effect. In the QED case, the light-by- 
light scattering graphs have an anomalous dependence on the external charges and a cut structure 
corresponding to particle production. In addition, we note that the non-analytic contributions to 
ay in higher orders in QCD arise from corrections to the "H-graph" . Therefore it could be argued 
that these types of diagrams should be excluded when defining the V-scheme in QCD and QED. 

4 We expect the main effects from including the quark masses at the one-loop level as this is 
the leading term in the ■(/'-function. However, at small scales the higher order terms will become 
important, especially since the relative importance of the Np term is larger for than for 
A study at the two-loop level requires the massive two-loop diagrams which is work in progress 
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FIG. 1. Single insertion of massive quark-antiquark loop into a gluon propagator, giving the 
quark part of the one-loop gluon vacuum polarization. 

n<°>(Q, Q , « ) = T fJ2^ (jf 6*(1 - z) In (1 + z(l - z) Pi {Q))- 

j\z(l-z) ln(l + z(l-z) Pi (Q ))] 

where p = Q 2 /m 2 , Tp = |, the superscript (0) indicates the one-loop order, the subscript q 
indicates the quark-part and the sum runs over all quarks (n). Thus the quark component 
of the one-loop ^-function is: 
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This gives^ the contribution to Np from quark flavor i, 

N F,v(pi) = 6 / t~ — n — r~ = 1 — + 3/2 ; tanh \ — rr> 8 

which is displayed in Fig. § as a function of p. Thus, by keeping the explicit quark mass 
dependence, N F becomes an analytic function of the scale Q. 
In fact, the approximate form: 

(l + £) 1 (9) 

gives an accurate approximation to the exact form to within a percent over the entire range 
of the momentum transfer^. 

The one-loop analytic Np t y is shown in Fig. [5] for various quark flavors (for reference, the 
quark masses (in GeV) we used are: m u = .004; m d = .008; m s = .200; m c = 1.5; m h = 4.5; 
m t = 175). 



5 This result was first obtained by Georgi and Politzer [39| in the MOM scheme and was applied 
to general gauge theories by Ross flicf . 

6 This approximate form can be obtained from using a rigorous double asymptotic series approach, 



knowing the behavior of the function at the low and high momentum transfer. |41] 
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FIG. 2. The curve shows the contribution to the continuous Npy for just one quark as a 

function of p = Q 2 /m 2 where m is the mass of the quark. Npy is found by using the massive quark 
part of the one-loop gluon propagator instead of using the theta function thresholds conventionally 
used in dimensional regularization schemes. 



We may now substitute the Npv into the one-loop QCD ip function coefficient: 



11 1 



2 3 

and thence into the QCD one-loop renormalization group equation for the coupling constant: 



dct\ 



dlnQ " , < 10 > 

We may then solve this renormalization group equation to yield an expression for ay which 
is analytic at mass thresholds. Note that the mass- dependence of the ip function applies 
specifically to the ay schemeQ. 



B. Commensurate scale relation between ay and a^g- 

We now relate the mass dependence of the ay scheme to the MS scheme using the 
commensurate scale relation [13|j36| between the two schemes. We use the NNLO results of 
Peter |33]. The first step is to invert Eq. (§) to obtain a-y^ as an expansion in ay, 



7 Given ay{Qo) one can obtain the coupling at other scales including the mass dependence by 
numerical iteration such as the fourth-order Runge-Kutta algorithm. 
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FIG. 3. Continuous JV>> for various quarks, lightest to heaviest goes top to bottom (d, c, b, t 
as one proceed downwards; the u and s plots are virtually identical at this scale to the d). Q runs 
from 1 to Mz GeV (for reference, the quark masses (in GeV) used are: m u = .004; = .008; 
m s = .200; m c = 1.5; m& = 4.5; m t = 175.) 



ans(Q) = MM) + mi , r (I) + m%v (|) sm. + . . . . 



The needed commensurate scale relation is obtained by fixing the scales M in Eq. ( |TT|) such 
that the and t/jW dependent parts of the coefficients miy and m2,y are absorbed into 
the running of the coupling ay(M). This insures that all vacuum polarization dependence 
is summed into the heavy quark potential. Application of this procedure in next-to-next-to 
leading order (NNLO), using the multi-scale approach P5fl , gives the following scale- fixed 
relation between ay and the conventional MS, 

am{Q) = a v {Q*) + lN c ai(Q ' 
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above or below the quark mass threshold^. The coefficients in the perturbation expansion 
have their conformal values, i.e. the same coefficients would occur even if the theory had 
been conformally invariant with ip^ = and thus do not contain the diverging (ip^ a s ) n n\ 
growth characteristic of an infrared renormalon The next-to leading order (NLO) 



coefficient \Nc is a feature of the non-Abelian couplings of QCD and is not present in QED. 
The commensurate scales Q* and Q** are given by 
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whereas to this order Q*** is not constrained. However, a first approximation is obtained by 
setting Q*** = Q**. Also note that Q* is unchanged when going from NLO to NNLO. The 
scale Q* arises because of the convention used in defining the modified minimal subtraction 
scheme. Comparing the scales Q and Q* we find that the scale in the MS scheme (Q) is a 
factor ~ 0.4 smaller than the physical scale (Q*). 

Alternatively, one can write the relation between a^s and ay as a single-scale commen- 

In this procedure Q* = Q** where 



surate scale relation 42 



Q* = Q exp 
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(15) 



The conformal coefficients are the same in the two procedures^. However, the single-scale 
form has the advantage that the non-Abelian perturbation theory matches in a simple way 
the corresponding Abelian perturbation theory in the limit — > with CfQ s arid Nf/Cf 
fixed 13]. For N c = 3 we have ln(Q*/Q) = 5/6 + A.178a v /n. 



C. Definition of the Analytic 

We now adopt the commensurate scale relation with the effective charge of the effective 
potential as a definition of the extended scheme ttM§"i 

MQ)-^) + ^^ + -> (is) 

O 7T 

for all scales Q. Eq. ( |T6"D not only provides an analytic extension of dimensionally regulated 
schemes, but it also ties down the renormalization scale to the physical masses of the quarks 



8 Note that the NNLO results depend crucially on whether or not the "H-graph" is included in 
the definition of the heavy quark potential since it is the unique source of the vr 4 A r ( 2. terms in the 
NNLO coefficient. We thank M. Peter for communications on this point. 

9 Both the multiple and single-scale setting methods generate a term proportional to CfNc in 
the NNLO conformal coefficient. The origin of this term, which has the same color factor as an 
iteration of the potential, is not clear and should be further investigated. 
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as they enter into the vacuum polarization contributions to ay. There is thus no scale 
ambiguity in perturbative expansions in ay or 5ms. 

Taking the logarithmic derivative of the commensurate scale relation given by Eq. (|l~6"|) 
with respect to In Q we can define the ^-function for the 5m§- scheme as follows, 



fe(Q) = MQ1 + 2 2N o caviQ ** ) MQn- 



7T 



(17) 



To lowest order this gives ^>M(Q) = ipy(Q*), which in turn gives the following relation 

(0) n „A aA°) 



between Npj^ and N FV , 
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where to lowest order, Q* = exp(5/6)Q- 

We can also use the approximate form given by Eq. 
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In other words the contribution from one flavor is ~ 0.5 when the scale Q equals the quark 
mass mj. Thus the standard procedure of matching aj^^i) at the quark masses is a zeroth 
order approximation to the continuous N F . 
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FIG. 4. The continuous -^vj^jg in the analytic extension of the MS scheme as a function of 
the physical scale Q. (For reference the continuous Np is also compared with the conventional 
procedure of taking A^^ to be a step-function at the quark-mass thresholds.) 

Adding all flavors together gives the total Np^(Q) which is shown in Fig. |4]. For reference 
the continuous Np is also compared with the conventional procedure of taking Np to be a 
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step-function at the quark-mass thresholds. The figure shows clearly that there are hardly 
any plateaus at all for the continuous Np^(Q) in between the quark masses. Thus there is 

really no scale below 1 TeV where Np^(Q) can be approximated by a constant. In other 
words, for all Q below 1 TeV there is always one quark with mass m, such that m 2 -C Q 2 or 
Q 2 3> mf is not true. We also note that if one would use any other scale than the BLM-scale 
for Np^—(Q), the result would be to increase the difference between the analytic Np and 
the standard procedure of using the step-function at the quark-mass thresholds. 



D. Comparing the Analytic a-^g{Q) with a^g- 

We can obtain the renormalization group equation for the analytic extension of the MS 
coupling 5ms by using NfL&(Q), etc.: 

da m (Q) _ r(o)/ n xa|s(Q) , om 

ding ~ tt + "'- {M) 

The solution to Eq. (p0|) provides an analytic scale- fixed extension of oc^s(fi), which we 
have denoted as 5m§-. The result can be compared with the standard method of computing 
q?ms, based on the evolution with distinct ip functions for different quark mass regimes. 
When doing this comparison one has to keep in mind that it is possible to take quark mass 
threshold effects into account also in the MS scheme when calculating an observable. In the 
next section we will compare our analytic extension of the MS scheme with the standard 
treatment of quark mass threshold effects for the hadronic width of the Z-boson. However, 
for most observables the quark mass threshold effects are not known and thus it is also 
important to compare an d Oms directly. 

Fig. H] shows the relative difference between the two different solutions of the 1-loop 
renormalization group equation, i.e. (ot^(Q) — cxus(Q))/ot^(Q). The solutions have been 
obtained numerically starting from the world average |2lJ a-Mg{Mz) = 0.118. The figure 
shows that taking the quark masses into account in the running leads to effects of the order 
of one percent, most especially pronounced near thresholds. 

In addition the figure shows the results obtained by using two different scales in Npl—(Q), 



namely 2Q and Q/2, when solving Eq. (20). This shows clearly that the BLM-scale min- 
imizes the difference between solutions using the continuous and discontinuous Np. The 
other two scale choices gives differences of several percent for small Q. 

We see from the figure that the effect of treating thresholds continuously can be of the 
order of a few percent in the magnitude of the QCD coupling when running down from Mz to 
m c . This is a significant difference, at the level of the precision of current a s determinations. 
The primary factor which influences the running is the value of the one-loop ip function, 
r = T — I-^f' a l ar g er value of Np gives a smaller ip^ and makes a s run more slowly; 
conversely a smaller value of Np gives a larger ip^ and a s runs more quickly. 

We can trace the difference in the couplings as follows: at Mz, the continuous function 
is above the discrete threshold value of 5, but goes below it at 30 GeV; it remains 
below the discrete threshold value until rrib where it becomes larger and remains larger until 
~ 3 GeV, where it becomes smaller again etc. 
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FIG. 5. The solid curve shows the relative difference between the solutions to the 1-loop renor- 
malization group equation using continuous Np, a^g(Q), and conventional discrete theta- function 
thresholds, a^g(Q). The dashed (dotted) curves shows the same quantity but using the scale 2Q 
(Q/2) in iVp-^g. The solutions have been obtained numerically starting from the world average 
1| a m (M z ) = 0.118. 

Thus, running down from Mz, «ms runs slower than a-^g until 30 GeV where the difference 
between them begins to close as 5ms runs faster than a^; at ~ 8 GeV the difference starts 
to increase again until the b quark threshold where starts to run slower than a^g and 
the difference between the two decreases until ~ 3 GeV, etc. this behavior forms the peaks 
seen m Fig. H Thus we see that aM§- will end up higher than when running down to low 
momentum transfers starting from Mz- 



IV. APPLICATIONS 

In this section we will show how to compute an observable using the analytic extension 
of the MS scheme and compare with the standard treatment of quark mass threshold effects 
in the MS scheme. The essential difference between the perturbative expansions in the a-M§ 
and 5ms couplings are terms that contain quark masses. In the analytic scheme the quark 
mass effects are automatically included whereas in the MS scheme they have to be included 
by hand for each observable. 

For some observables, such as the hadronic width of the Z-boson and the r lepton semi- 
hadronic decay rate, corrections due to non-zero quark masses have been calculated within 
the MS scheme [p!8|-|20|jm . To be specific we are interested in the so called double bubble di- 



agrams where the outer quark loop which couples to the weak current is considered massless 
and the inner quark loop is massive. Other types of mass corrections, such as the double 
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triangle graphs where the external current is electroweak, are not taken into account by the 
analytic extension of the MS scheme. (For a recent review of higher order corrections to the 
Z-boson width see 

To illustrate how to compute an observable using the analytic extension of the MS scheme 
and compare with the standard treatment in the MS scheme we consider the QCD corrections 
to the quark part of the non-singlet hadronic width of the Z-boson, T^ d q . Writing the QCD 
corrections in terms of an effective charge we have 



-pTVS 
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(21) 



where the effective charge a N ^{s) contains all QCD corrections, 
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(22) 



The functions F and G are the effects of non-zero quark masses for light and heavy quarks, 
respectively. In the following we will not restrict ourselves to the case yfs = Mz since we 
want to compare the two treatments of masses for arbitrary s. Thereby the number of light 
flavors Nl will also vary with s. We will also assume that the matching of Np is done at 
the quark masses. Thus a quark with mass m < fi is considered as light whereas a quark 
with mass m > fi is considered as heavy. 

To calculate a^(s) in the analytic extension of the MS scheme one first has to apply 
the BLM scale-setting procedure which absorbs all the massless effects of non-zero Np into 
the running of the coupling. This gives, 
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(23) 



where 



Q* = exp 



,11 2, 



^ = 0.7076x/s. 



(24) 



Operationally, one next simply drops all the mass dependent terms in the above expression 
and replaces the fixed Np coupling atj^ with the analytic 5m§". (For an observable calculated 
with massless quarks this step reduces to replacing the coupling.) In this way both the 
massless Np contribution as well as the mass dependent contributions from double bubble 
diagrams are absorbed into the coupling and we are left with a very simple expression, 



a 



NS, 



s) 



71 



71 



(25) 



This simple expression reflects the fact that the effects of quarks in the perturbative coeffi- 
cients, both massless and massive, should be absorbed into the running of the coupling. 
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To compare with the ordinary MS treatment we need the functions F and G in Eq. ( |23| ) . 
Expansions in terms of m 2 /s and s/m 2 can be found in |TB|JTP| ,P[] whereas they have been 
calculated numerically in [20]. In addition the a% correction due to heavy quarks has been 
calculated as an expansion in s/m 2 in M. It should also be noted that the function G was 



first calculated for QED |4"5 |. Here we will use the following expansions, 
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(26) 
(27) 



which are good to within a few percent for m 2 /s < 0.25 and s/m 2 < 4 respectively. We will 



also use the relation |20|, 





6 \ s V 12 3 s 



(28) 



to get F in the interval 0.25 < m 2 / s < 1 since the expansion of F in terms of m 2 /s breaks 
down for m 2 /s > 0.25. 

Before carrying out the comparison of the analytic extension of the MS scheme with the 
standard treatment it is instructive to look at the effective contribution to a^(s) from one 
flavor with mass m as a function of s. To make the arguments more transparent we will use 
the renormalization scale fi = yfs when doing this. For small s, when the quark is considered 
heavy the contribution is given by G(m 2 / s) whereas for larger s the quark is considered as 
light and contributes with F{m 2 /s) — t| + -£3. Normalizing to the massless contribution 



11 
12 



+ 1^3 gives the contribution to the effective Np in the ag -coefficient, 
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for a/s > m 



(29) 



which is shown in Fig. || as a function of s/m 2 . 

At first it might seem unnatural that the effective contribution to Np in the ttg-coefficient 
is negative for heavy quarks. However, this is a characteristic feature of the standard MS 
scheme which arises from the fact that the number of flavors in the running of the coupling 
is kept constant. Starting from a scale well below the threshold the number of flavors in the 
running as well as in the a 2 coefficient is not affected by the heavy quark. As the threshold 
is approached from below the number of flavors in the running should increase which would 
make the running of the coupling slower (since ip^ would be smaller) which in turn should 
lead to a larger ckms- But, since the number of flavors is kept constant in the running this 



15 



N F e 'rs(s/ m 2 ) 



0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 

10~ 2 10~ 1 1 10 10 2 

s/m 2 

FIG. 6. The effective contribution to Np in the aig-coefficient in the standard MS scheme from 
a quark with mass m as a function of s/m 2 (using fi = y/s). The discontinuity between the two 
expansions in s/m 2 and m 2 / s can be seen at the nonanalytic point s/m 2 = 4. 

effect has to be taken into account by adding a positive contribution to the ag-coemcient, i.e. 
the function G. Since the massless contribution is negative this means that the contribution 
to Np becomes negative for a heavy quark. Once the threshold has been crossed the number 
of flavors in the running changes and the need to compensate for a too small vanishes 
rapidly as the scale is increased above the threshold. For scales well above the threshold 
the mass-effects are negligible and the massless result is regained as F goes to zero. This 
should be compared with the analytic MS scheme where Np is increased continuously in the 
running. 

To compare the analytic extension of the MS scheme with the standard MS result for 
a r,q( s ) we wn l a PPly t ne BLM scale-setting procedure also for the standard MS scheme. 
This is to ensure that any differences are due to the different ways of treating quark masses 
and not due to the scale choice. In other words we want to compare Eqs. (|23"D and (|25|) . As 
the normalization point we use a^(Mz) = 0.118 which we evolve down to Q* = 0.7076M^ 
using leading order massless evolution with Np = 5. This value is then used to calculate 
a^(Mz) = 0.1243 in the MS scheme using Eq. (H). Finally, Eq. (H) gives the normaliza- 
tion point for <5ms((3*). 

Fig. [7| shows the relative difference between the two expressions for a^g(s) given by 
Eqs. (^) and fl25|) respectively. As can be seen from the figure the relative difference is 
smaller than 0.2% for scales above 1 GeV. Thus the analytic extension of the MS scheme 
takes the mass corrections into account in a very simple way without having to include an 
infinite series of higher dimension operators or doing complicated multi-loop diagrams with 
explicit masses. 
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FIG. 7. The relative difference between the calculation of apf ( s ) ™ the analytic extension of 
the MS scheme and the standard treatment of masses in the MS scheme. The discontinuities are 
due to the mismatch between the s/m 2 and m 2 /s expansions of the functions F and G. 

V. CONCLUSION 

An essential feature of the av(Q) scheme is that there is no renormalization scale am- 
biguity, since Q 2 is the physical momentum transfer. The ay scheme naturally takes into 
account quark mass thresholds, which is of particular phenomenological importance to QCD 
applications in the intervening mass region between those thresholds. In this paper we have 
utilized commensurate scale relations to provide an analytic extension of the conventional 
MS scheme in which many of the advantages of the ay scheme are inherited by the a-^g 
scheme, but only minimal changes have to be made to the standard a-^g scheme. Given 
the commensurate scale relation, Eq. fljjD , connecting a-^ to ay expansions in are 
effectively expansions in ay to the given order in perturbation theory at a corresponding 
commensurate scale. Unlike the conventional scheme, the modified 5m§- scheme is an- 
alytic at quark mass thresholds, and it thus provides a natural expansion parameter for 
perturbative representations of observables. In the Abelian limit N c — > 0, 5^ scheme 
agrees with the standard effective charge method of QED. 

We have found that taking finite quark mass effects into account analytically in the 
running, rather than using a fixed Np between thresholds, leads to effects of the order of 
one percent for the one-loop running coupling, with the largest differences occurring near 
thresholds. These differences are important for observables that are calculated neglecting 
quark masses and could in principle turn out to be significant in comparing low and high 
energy measurements of the strong coupling. 

We have also found that our extension of the MS scheme, including quark mass effects 
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analytically, reproduces the standard treatment of quark masses in the MS scheme to within 
a fraction of a percent. The standard treatment amounts to either calculating multi-loop 
diagrams with explicit quark masses or adding higher dimension operators to the effective 
Lagrangian. These corrections can be viewed as compensating for the fact that the number 
of flavors in the running is kept constant between mass thresholds. By utilizing the BLM 
scale setting procedure, based on the massless Np contribution, the analytic extension of the 
MS scheme correctly absorbs both massless and mass dependent quark contributions from 
QCD diagrams, such as the double bubble diagram, into the running of the coupling. This 
gives the opportunity to convert a calculation made in the MS scheme with massless quarks 
into an expression which includes quark mass corrections from QCD diagrams by using the 
BLM scale and replacing a-^ with 5ms- 

For simplicity we have analyzed the mass corrections arising from analyticity only to 
leading order in QCD. For further precision, our analysis will need to be systematically 
improved. For example, at higher orders the commensurate scale relation connecting ay to 
ojms will have to be corrected with finite mass effects. We have seen that the BLM scale 
minimizes the difference between the analytic and the conventional aMg-coupling. Thus, 
these kinds of corrections are not likely to decrease the difference between the analytic and 
the conventional ^Ms-coupling. 

Finally, we note the potential importance of utilizing the ay effective charge or the 
equivalent analytic 5ms scheme in supersymmetric and grand unified theories, particularly 
since the unification of couplings and masses would be expected to occur in terms of physical 
quantities rather than parameters defined by theoretical convention. 
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